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. We show that the formal procedure of integrating out the degrees of freedom of the small spheres 

^— ' in a binary hard-sphere mixture works equally well for non-additive as it does for additive mix- 

tures. For highly asymmetric mixtures (small size ratios) the resulting effective Hamiltonian of 
^-^i the one-component fluid of big spheres, which consists of an infinite number of many-body inter- 

im ' actions, should be accurately approximated by truncating after the term describing the effective 

pair interaction. Using a density functional treatment developed originally for additive hard-sphere 
mixtures we determine the zero, one, and two-body contribution to the effective Hamiltonian. We 
demonstrate that even small degrees of positive or negative non-additivity have significant effect on 
the shape of the depletion potential. The second virial coefficient B2, corresponding to the effective 
pair interaction between two big spheres, is found to be a sensitive measure of the effects of non- 
y j — ^ , additivity. The variation of B2 with the density of the small spheres shows significantly different 

Q ■ behavior for additive, slightly positive and slightly negative non-additive mixtures. We discuss the 

' possible repercussions of these results for the phase behavior of binary hard-sphere mixtures and 

suggest that measurements of B2 might provide a means of determining the degree of non-additivity 
in real colloidal mixtures. 
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Mixtures of hard spheres play a pivotal role in the statistical mechanics of liquids. Not only do they provide a 
realistic reference system for describing the structure and thermodynamics of simple atomic mixtures and mixture 
of colloidal particles, they are also of considerable intrinsic interest. In particular, investigating the properties of 
asymmetric binary hard-sphere mixtures became a topic of much activity when it was recognized by Biben and 
Hansen that such athermal mixtures might afford important examples of pure entropy-driven fluid-fluid phase 
separation. The most studied model is that of additive hard spheres, where the cross diameter <ti,s = (cTbb + ass)/'2 
and cTbb refers to the big-big and (Tss to the small-small diameters. There is now compelling evidence to suggest that 
additive mixtures, with sufficiently small size ratios ass/o'tb-, do undergo fluid-fluid phase separation but this transition 
remains metastable w.r.t. the fluid-solid transition The other well studied model is the so-called Asakura-Oosawa 
(AO) model of a colloid-polymer mixture in which the colloid-colloid interaction is hard-sphere like, with diameter 
Ch ' cTcc, and the colloid-polymer interaction is also hard with diameter Ccp, whereas the polymer-polymer interaction is 
, zero, i.e. app — 0, corresponding to ideal interpenetrating coils. The cross diameter dcp — {cFcc + '^Rg)l'^, where Rg is 
(— I ' the radius of gyration of the polymer. Various approximate theories [p|^ and some simulation studies for simplified 
Q , versions of the AO model showed that when the size ratio 2Rg/c7cc is larger than about 0.35 fluid-fluid phase 
O ■ separation is stable w.r.t. the fluid-solid transition. The AO model can be regarded as an extreme non-additive 
^ ' mixture with acp > (cTcc + o'pp)/2. 

One can, of course, investigate mixtures with arbitrary non-additivity and recent studies of binary hard-sphere 
mixtures have indicated that a small degree of positive non-additivity in the cross diameter ats might be sufficient to 
induce a fluid-fluid transition P,p^ . By employing an effective one-component treatment and a variety of liquid state 
perturbation theories, Louis et. al pT] have demonstrated that non-additivity should have a profound effect on both 
the fluid-fluid and fluid-solid transition in highly asymmetric hard-sphere mixtures. However, their treatment is based 
on an empirical approximation [ p^ for the effective (depletion) potential between two big spheres rather than any 
systematic derivation of an effective one-component Hamiltonian for the non-additive hard-sphere mixture. In this 
paper we develop such an approach following the path that was trodden in the recent studies of the additive mixture 
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1^ and in the special case of the AO model . We show that the same formal technique of integrating out the degrees 
of freedom of the smaller species, used in Ref. 1^, applies equally well to non-additive mixtures (see Sec. II A ). 



We were motivated toward such an approach by the following considerations: a) treating highly asymmetric mixtures 
by brute force simulation is beset by ergodicity problems and slow equilibration when the packing fraction of the small 
species is substantial, b) for small size ratios three and higher body potentials in the effective Hamiltonian do not 
have a significant effect on the phase behavior of the additive or of the AO [Qjs) model, i.e. in both models phase 
transitions are determined primarily by the effective pairwise potential between the big particles and we expect the 
same to be true for small size ratios in systems with intermediate degrees of non-additivity. c) the effective pairwise 
potential which arises in the formal development of the theory is the depletion potential, introduced into colloid 
science by Asakura and Oosawa, and now much studied by theory, simulation |]l3|,|l4| and experiment. Indeed a 



variety of experimental techniques |15| have been developed to measure the depletion potential between a colloidal 
particle, immersed in a sea of small colloids or non-adsorbing polymer, and a wall or another big colloid. Interpreting 
the results of such experiments requires a reliable theory. Recently we have shown that a density functional 

approach (DFT) provides an accurate means of calculating the depletion potential for additive hard spheres. In the 



present paper we show (Sec. [IB ) that the same DFT approach remains valid for the non-additive case and can 
therefore be used to investigate a much wider class of depletion potentials than one might have suspected a priori. We 
find that even very small degrees of non-additivity can have a very pronounced effect on the shape of the depletion 
potential, which leads to significant changes of the magnitude, and sometimes the sign, of the second virial coefficient 
B2 associated with the total effective potential between two big spheres. Such changes in B2 may, in turn, have 
repercussions for the phase behavior, and may be directly accessible by experiments. 

We begin by defining what we mean by non-additivity. A non-additive binary hard-sphere mixture is characterized 
by the diameters (distances of closest approach of the centers of the spheres) ai,{,, ass, and ai,s where the subscripts b 
and s denote big and small. These diameters describe the (pairwise) interaction potentials between two spheres 



$ fr) = <^ ^^'1 (1) 



00 r < ai 
"ijy J — Q otherwise, 

with i,j £ {s, b}. We follow the usual convention and introduce the non-additivity parameter A via 

abs = ^{<Jbb + crss){l + A) (2) 

and allow A to be positive or negative. Additive hard-sphere mixtures have A = 0. 

In the following we consider two different routes to introducing non-additivity into a binary hard-sphere mixtures. 
In both routes we keep the diameter of the big spheres ahb and the number density of the small spheres in the 
reservoir fixed but allow for (a) changes in the diameter of the small spheres ass while keeping the cross diameter 
CTbs constant or (b) changes in the cross diameter ct^s while keeping the diameter of the small spheres ass constant. 
These two types of changes are illustrated in Fig's |^ and ^ respectively. 

In the first route (a) , which was introduced in an earlier Letter by two of us Jlq] , we can smoothly follow a path 
that connects the additive hard-sphere mixture to the Asakura-Oosawa model |^^. As mentioned above, the small 
particles are modeled as ideal gas particles with zero diameter aff = and A = A^'^ > 0, keeping the same cross 
diameter as in the additive case, i.e. we require 

<7bs = liTbb + cTss{A)){l + A) = const. = <f = ^{ahb + <f ) (3) 

SO that 

—add — A 

^..(A)- "-^;7^ . (4) 

The AO limit, defined by (Tss(A"^'^) — 0, implies A"^*^ = a'^s'^/ahb = q, the (fixed) size ratio. Within route (a) 
choosing A > A^'^ would give rise to iTss(A) < which is unphysical, of course. For a given value of the non- 
additivity parameter A and a fixed number density of the small spheres Ps it follows that their packing fraction 
r?^(A) — 7rcr^^(A)p^/6 also varies with A: 

addf 9 - A 



where 77°'^'^ = ^/^(O). Clearly ?yJ(A) decreases with increasing A in the range < A < A'^*^ = q. On the other 
hand, for negative non-additivity, A < 0, r]l{A) increases rapidly with increasing |A| and for studies of the fluid 
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phase we should restrict ??^(A) < rj^^'^'^^ = 0.494, the value of the packing fraction at the bulk freezing transition of 
hard-spheres. 

If we consider route (b), the route used by most previous authors both ai,i, and ass are kept fixed and 

the cross diameter ats becomes a function of A through Eq. (||). Note that this case is very different from case 
(a) as the packing fraction rj^ of the small spheres remains constant through the whole range of non-additivity and 
the Asakura-Oosawa limit can only be reached in the limit of ^ in which the depth of the depletion potential 
approaches zero as well. 

Examples of depletion potential and second virial coefficients calculated from both routes are given in Sec. HI. We 
conclude in Sec. IV with a summary of our results and a discussion of their possible relevance for experiments. 



II. MAPPING A BINARY MIXTURE ONTO AN EFFECTIVE ONE-COMPONENT SYSTEM 

In this section we map the binary mixture of non-additive hard spheres onto an effective one-component system 
of big particles. To this end we begin by formally integrating out the degrees of freedom of the small spheres and 
determine the form of the effective Hamiltonian of the one-component fluid. We demonstrate that each term of the 
effective Hamiltonian of the non-additive system can be determined using a theory for an additive mixture. In a 
second step we calculate explicitly the leading contributions to the effective Hamiltonian for the binary mixture of 
non-additive hard spheres using a DFT designed for the additive hard-sphere mixture. 



A. Formal mapping for additive and non-additive mixtures 

We follow the procedure developed in Ref. |^ to formally integrate out the degrees of freedom of the small spheres 
in a homogeneous mixture of Ni, big and Ns small spheres in a macroscopic volume V. The Hamiltonian of our 
mixture is given by 

H = K + Hbb + Hss+Hbs, (6) 

with K the total kinetic energy of the mixture, leading to a trivial contribution to the free energy, and three potential 
energy contributions: 

Nb N, JVi, N, 

i<j i<3 i=l i=l 

where the pairwise hard-sphere interaction potentials $ij are defined in Eq.(|l|). y\, denote the coordinates of big 
particle i and small particle respectively. It was shown in Ref. that the thermodynamic potential F(Nb, Zs, V) 
of a general binary mixture in the semi-grandcanonical ensemble can be written in terms of an effective Hamiltonian 
H'^-f^ via the relation 

exp(-/3F) = Tnexpi-PH^ff), (8) 

where /3 = l/fc^T, and Zs = Aj'^ exp(/3/is) is the fugacity of the (small) species s, fixed by the reservoir. Tr;, is 
shorthand for I/N^IA^,^" times the volume integral over the coordinates of species v. In Eq. (||) it denotes the 
classical trace over the degrees of freedom of the (big) species b. The effective Hamiltonian H"^^^ is given ||] by 

H'^^f = Hbb + n, (9) 

with ^l = ^l{Nb, Zs, V; {r''}), the grandcanonical potential of the fiuid of species s subjected to the external potential 
of a fixed configuration {r''} of Nb particles of species b, defined by 

oo 

exp(-/3r!) = exp{PtisNs)TTs exp(-/3(F,, + Hbs)). (10) 

N,=0 

Using a Mayer cluster expansion it was shown ||^ that can be written as a sum of terms r2„ which describe the 
simultaneous interaction of n particles of species b with the 'sea' of species s, i.e. 
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n=0 



(11) 



This result is valid for arbitrary (integrable) pairwise potentials. We emphasize that in Eq. (||) all direct interactions 
between species b, the big particles, are contained in Hi,b while Q describes interactions between species s, small 
spheres, and between the big and small ones. It is precisely this separation into a term which contains just the big-big 
interactions and those which contain small-small and big-small interactions which allows us to calculate the leading 
terms r2„ for non-additive hard-sphere mixtures using a theory for an additive mixture. More specifically, for a given 
fixed configuration of large particles, is completely determined by the parameters Zs, ctss and abs- In other words 
all the terms fin would have an identical form for an additive or a non-additive system. The only difi'erences arise in 
Hbb, where abb constrains the possible positions of the big particles to |rj — r^| > abb for all i,j. For the additive case 
abb = '2abs — Ctss, whilc in the non- additive case abb can vary more widely depending on the value of A. 



1. Zero body term Q.q 

The first term, i7o, is the grandcanonical potential of a sea of small spheres with fugacity Zs without any big sphere 
present and it follows that 

no{zs,V) ^ ~Ps{zs)V, (12) 

with Ps{zs) the pressure of the reservoir of small spheres. Since this term is intrinsic to the small-sphere fluid it is not 
affected by introducing non-additivity. 



2. One body term fii 

For a homogeneous system the one-body term f2i is of the form 

ni{Nb,Zs)=NbLOi{zs) (13) 

with |] 

exp(-/3c.i(z,)) = (exp(-/3H,W)),^ (14) 

where h'^^ denotes the interaction between Ns small spheres and n > \ big spheres and the brackets (. . .)z^ refer 
to an ensemble average in the reservoir of small spheres, wi(zs) can be identified as the difference in grandcanonical 
potential between a sea of small spheres at fugacity Zs with and without a single big sphere. By considering the 
otential distribution theorem and the definition of the one-body direct correlation function c\^^ (r**) one can show 
"9|,p^ that this difference in grandcanonical potential can be expressed as 

/3wi(z,) = - lim 4'^(oo), (15) 

where c"fP {oo) denotes the direct correlation function of a big sphere evaluated in the bulk mixture. The limit 
/if, — *■ — oo implies that the chemical potential ^b of species b is made sufficiently negative that only one big sphere 
is present. c[^^(oo) is proportional to the excess (over ideal) chemical potential of species &, i.e. — c[^^(oo) = (3^1^ 

and, in general, depends on the density of both species. However, in the limit /if, — )■ — cx), c\p {oo) depends only on 
the fugacity Zg. Since only one big sphere is involved [this is explicit in Eq. (p^)] non-additivity plays no role in 
determining 17 1. One merely specifies ass and then abs describes the interaction between small spheres and a fixed 
big one. 

As noticed in Re f. |po| , the one-body term uji{zs) determines the Henry's law constant }i{zs) of the fluid. The latter 
can be defined by pl[: 

h{zs) = lim (16) 

Pb^O Zb{pb) 



where Zb — h^^^ exp(/3/if,) is the fugacity of species b. It follows that 
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h{zs) = exp(-/3wi(z^)). (17) 

h(zs) does not depend on b-b interactions; deviations of h{zs) from unity reflect the average effect of b-s interactions 
at a fixed fugacity Zg of small spheres. Thus, given some means of calculating cjj^\oo), in the limit pb 0, one can 
determine the Henry's law constant. 



3. Two body term Q2 



The two body term Q2 is given by: 



n.iNb, z,; {r"}) = 5]u;2(^.; |r? - r% (18) 



i<j 



where the pair potential uj2 is defined by j2j 



exp(-/3c.,(z.;rM = ^'"^^ ^^f'!^?" (19) 



(exp(-/3<)))2^ 



with = |r^ — r^| the distance between the centers of the big spheres. Equivalently we can use the variable 

h = r\j — <Jbb^ i.e. the distance between the surfaces of two big spheres. uj2{zs;r\j) = W{h — r'lj — atb) is the grand 
potential difference between a sea of small spheres, at fugacity Zg, containing two big spheres at finite separation 
and one in which the separation r^^- = 00. Although the positions of two big particles are involved, the big-big 
interaction does not enter explicitly into the calculation of uj2, the diameter abb merely acts as an external parameter 
which restricts the minimum separation to r^j — abb- In other words, for a given Zg, abs and ass: the calculation of 
L1J2 is the same for both additive and non-additive systems. 

UJ2 can be identified with the well-known depletion potential between the two big spheres PJl6| and expressed in 
terms of the one-body direct correlation function: 



Mzs-.v) = hm (cW(oo)-cW(r)) (20) 

neous situation in whi 
small spheres and a big (test) particle is inserted at r fl^ . 



where c^^^(r) refers to an inhomogeneous situation in which a big sphere fixed at the origin exerts its field on the 



4. Three body and higher order terms Q,n>3 



The three-body term can be written as a sum of three-body potentials 0^3(23; r^ ^,) which can, in turn, be 



expressed in terms of ensemble averages (exp(— /^i?^^"''))^^ with n ~ 1,2,3 ||^. Once again big-big interactions are 
not involved in the calculation and abb simply specifies the physically allowed configurations of the three big spheres. 
Clearly non-additivity plays no role. The same argument applies for the higher-body [n > 3) contributions to Q, 
although the specification of the allowed configurations of the big spheres becomes increasingly complicated as n 
increases. 

In practice the calculation of Qn for n > 3 is tedious and determining the phase behavior for effective Hamiltonians 
which include these and higher-body interactions would be very cumbersome. Contributions to the many body terms 
arise from two mechanisms: (1) Directly, if for a given abb, <^bs is large enough to allow the overlap of more than 
two depletion layers. (2) Indirectly, if correlations between the small particles, present for all non-zero ass, induce 
interactions between more than two particles. Contributions from mechanism (1) to r2„ with n > 3 are identically 
zero for {2abs — <Jbb) / <7bb < 2/-\/3 — 1 = 0.1547 0, while for {2abs — o'bb)/<7bb < •\/3/2 — 1 = 0.2247 the contributions 
to fin with n > 4 terms are zero p2| etc... If ass > 0, then mechanism (2) will induce additional contributions to the 
many-body terms at all values of abs and abb- 

For the extreme non-additive AO model, where ass = but abs is non-zero, only mechanism (1) contributes to the 
many-body terms. An explicit form for the three-body term can be calculated |^3|, but this is still very tedious to 
evaluate. 

For additive binary hard-sphere mixtures with q = ass/o'bb < 0.1, where only the indirect mechanism (2) contributes, 



three-body contributions seem to be small ||l3|] . Recent DFT calculations 1 23 2j| of the three-body potentials show 
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that these are still much smaller than the two-body potentials for q = 0.2, where both mechanisms contribute, and 
that the indirect mechanism (2) can have an important effect on their shape. 

There is strong evidence from direct simulation studies of the additive binary system that retaining only two-body 
contributions in the effective Hamiltonian provides a very good account of the equilibrium phase behavior for q < 0.2 

. Similar conclusions were reached for a lattice version of the AO model [|j . In this study we shall focus on such 
highly asymmetric systems and neglect three and higher-body contributions to the effective Hamiltonian. With this 
assumption it follows that the structure of the homogeneous fluid (equilibrium correlation functions of the big spheres) 
is determined solely by the effective pairwise potential 

^eff{zs;r) ^<Pbb{r)+uj2{zs;r), (21) 

since fto and ili do not depend on the coordinates of the big particles ||2^. It is also straightforward to show that the 
phase equilibria of the binary mixture do not depend on the zero and one-body term . However these two terms do 
influence the total pressure and compressibility of the mixture . 

To complete this sub-section it is necessary to explain how to convert from p^, the number density of the small 
spheres in the reservoir to ps, the actual value of the small sphere density in the mixture. The average number of 
small spheres in the mixture is given by the thermodynamic relation 

dfis \ dlnzs I 

In Ref. iQl it was shown that an accurate approximation for the conversion can be obtained for additive hard-spheres 
with high asymmetry, q < 0.1, by truncating the expansion after the one-body term so that the average number of 
small spheres can be evaluated approximately from the formula 



\ a In 



The required density ps — {Ns)z^/V. We shall re-visit this approximation in the next subsection. 



B. Evaluation of flo, ^i, and ^2 within DFT 

1. Rosenf eld's fundamental measure DFT 



In order to make the mapping presented in Sec. II A explicit for the model of interest, namely the binary mixture of 
non-additive hard spheres, we calculate the terms fio, S^i, and U,2 within the framework of Rosenfeld's fundamental 
measure DFT p5[ | - a theory constructed for additive hard-sphere mixtures. As mentioned earlier, the reason why 
we can apply a theory constructed for additive binary mixtures is that only the interaction potentials ^ss, between 
two small spheres, and $f,s, between a big and a small sphere, enter into the calculation of ri„. This argument 
substantiates further the intuitive picture presented in Ref. fl^ . 

What is required for calculating fJi and ^2 is some means of determining c"^^ in the limit of vanishing density of 
big spheres. DFT provides a suitable route since 



4^^(r)^- /^-^y^ (24) 



where J^ex[pbT Ps\ is the excess (over ideal) intrinsic Helmholtz free energy functional of the mixture | [i6[ |. Thus given 
some prescription for the mixture functional one can calculate all the necessary ingredients. Rosenfeld's fundamental 
measure theory p5[ supplies an approximate functional T^x for an additive mixture of hard spheres of the form 



PTex[pb{v),p,{v)]^ J d'r^iM), (25) 

with weighted densities which depend on the fundamental geometrical measures of the spheres constituting the 
mixture. There are four scalar and two vector weight functions w^, with 1 < a < 6. Details can be found in Ref. [ p5| ; 
the weights depend on the radii Ri of each species. For a binary mixture the weighted densities are 
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n„(r)= ^ /dVp,(r')wL(r-r'). (26) 

i—s^b 



It is imp ortant to realize that once the reduced free energy density ^{{ua}) is specified the mapping described in 



Sec. [I A is completely determined within this (approximate) DFT framework. We choose to apply the original 



Rosenfeld functional [E5| 

* {"a! = -"oln(l - na) + + ^TTTi ^2-■ 27) 

(2) 

Recall that the two-body direct correlation functions c)^ , with i, j e b, s, obtained by taking two functional derivatives 

of this functional reduce to the Pcrcus-Yevick (|r — r'|) for a homogeneous hard-sphere mixture. The Rosenfeld 
functional has proven to be extremely successful in describing the structure of the inhomogcncous fluid phases of 
hard-sphere mixtures. If solid phases are to be considered, modifications to the original Rosenfeld functional should 
be made [ETI. 



S. Calculating Slo from DFT 

We begin the explicit mapping by noting that the equation of state underlying the Rosenfeld functional ||2^] is 
Ppy, the Percus-Yevick compressibility equation of state for additive hard-sphere mixtures. In order to calculate the 
zero-body term, however, we need only the equation of state for a one-component fluid so we set 



(1 - ?7P^ 



with ril = 7r(T^^p^/6, the reservoir packing fraction of the small spheres, flo follows directly from Eq. 



3. Calculating Q.i from DFT 
Using the Rosenfeld functional the direct one-body correlation function c^^'' can be written as 

cl^\^)--Y.Jd'^'{^^^)^,<i^-^')- (29) 

In the dilute limit in which the density of the big spheres — > 0, the weighted densities depend only on the density 
profile Ps(r) of the small spheres: 



n 



d./«te(j.) ^ / ^3^, p,(r')w^(r - r'). (30) 



Note that ps (r) then corresponds to a one-component fluid of species s . In bulk the density profile of the small 
spheres Ps(r) is constant and equal to p^, the vector weighted densities fta vanish and scalar weight functions reduce 

to those of a one-component fluid, i.e. 713 rj^, n2 — > Grjl/ass, ni Sr/^ / (Trcr^g ) and uq Qrfg/ {TTa^^). c[^''(oo) can 
be evaluated explicitly and the result expressed as 

/3c.i(z,)= hm -c«(cx3) = ^/3p,(z,)+47ri?2/?^(z,)-f :^-^-ln(l-r?:), (31) 

where the first term corresponds to the partial derivative d/dn^, the second to d/dn2 and so on. As the first term is 
proportional to the pressure of the reservoir of small spheres, given by Eq. ( [28| ) and the second term is proportional 
to the planar 'surface tension' 

. ,i|ii±i£> , (32, 

these have a natural interpretation as PsAT^ and ^AA terms, respectively. Although it is more difficult to give a 
physical interpretation of the last two terms in Eq. (|3T|) they are important, for example in obtaining the correct low 
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density limit - see below. Since -/3-i4^^(oo) is the excess chemical potential of species 6 in a uniform (bulk) mixture 
we can also determine this quantity starting from the bulk excess free energy density, differentiating w.r.t. the bulk 
density and then taking the limit ~> 0. The result is identical to that in Eq. (|3l|). 

Although we have derived this result starting from a theory developed for an additive mixture once we have a taken 
the limit pb — *■ the big-big interaction is not relevant. It follows that for an arbitrary non-additive mixture, Rh 
entering Eq. (31) should be defined as Rh = Ubs — crss/2. In the particular case of an additive mixture Rb reduces to 

CTbb/2. 

Note that the form of Eq. ( pT| ) is the same as that given by Henderson [|l9| in a scaled particle theory for cj^\oo). 

However, in Henderson's treatment c^^^ (oo) is expressed in powers of the variable R = ai,s- If one converts his Eq. (54) 
to an expression in terms of Rb = R—Ossl'^ one recovers •precisely Eq. (^l]). In other words, his scaled particle analysis 
agrees completely with the present DFT approach, attesting further to the consistency of the latter. Note also that in 
the low density Hmit — > 0, Eq. (|3l|) imphes /?wi(zs) —* 47rcr^^Pg/3 which is the exact limiting value - see Appendix. 



^. Converting from the reservoir packing fraction rjl to the packing fraction ps in the mixture 



Given explicit formulae for VLq and VLi we can employ Eq. ( p3[ ) to obtain an explicit conversion between r]s, the 
packing fraction of the small spheres in the system, and 77^, the packing fraction in the reservoir. To this end we first 
determine the fugacity Zg = exp(/3/ig^) of the small spheres within the framework of the Rosenfeld DFT approach. 
Pg^, the excess chemical potential of the pure small sphere fluid, is given by the Percus-Yevick (compressibility) result 
so that 

/i4,r-i3(,r)- + 5(,;)3y ^^^^ 



Using the expressions for JIq and ili given in this section together with Eqs. (23) and ( |33D we obtain for the packing 
fraction of small spheres in the system ris = Trag^ps/G. 

where q^ff = ags/'^.Rb is the effective size ratio. In the limit 77^ — > this result reduces to 

llsh: = {l~Pb{l + qefff) (35) 

which is the standard excluded volume expression, appropriate to an ideal gas of small particles However, for 
non-zero densities of small spheres Eq. ( |3^ ) predicts a non-linear dependence of 77s on p^.. This is illustrated in Fig. ^ 
where we plot ps versus 77^ for additive hard-sphere mixtures with (a) q = 0.1 and (b) q = 0.05. Our results are 
compared with those of direct Monte-Carlo simulations of the binary mixture j^]. The agreement between theory and 
simulation is excellent, with a small deviation occurring at pb = 0.74 and 0.10 for q = 0.1. Note that the free-volume 
theory of Lekkerkerker and Stroobants [2^ , which asserts that Ps/pl = (^{zs = 0;?7b), where a{zs = 0',Pb) is the free 
volume fraction evaluated for zero fugacity of the small spheres, significantly underestimates Ps/Ps at higher values 
of 77J We conclude that retaining only the two leading terms flo and Sli and employing PY theory for these 
quantities provides an accurate approximation for the free volume fraction 773/77^, at least for small values of q. 

For completeness we should mention that the theory employed in Rcf. g| to calculate Ps/Ps used the Carnahan- 
Starling result for Zs and Henderson's expression for cl^\oo) but with an empirical modification of the leading (Rl) 
term. It is now clear that there was no need to make such a modification; the confusion arose from the improper 
identification of the parameter R in Henderson's theory. Fortunately the numerical results presented in Fig. 13 of 
Ref. 1^ are very close to those given by the present, fully consistent theory. For future applications we recommend 
that Eq. (ph should be used. 



5. Calculating Q2 from DFT 

The two-body contribution ^2, given in Eq. (|l8|), requires the calculation of the depletion potential U!2{zs; r) given 
by Eq. (|2^) . This can be carried out using the procedure described in Ref. |l^] . We first calculate the (inhomogeneous) 
equilibrium density profile Ps{r) of the small spheres near a fixed big sphere of radius Rb- This is used in Eq. (|3^) to 

determine the relevant weighted densities which then determine c^\r) via Eq. (|29|). 

We emphasize that the mapping of a depletion potential in a non-additive system onto one in an additive mixture 
is exact and has been applied in a recent study of generalized effective potentials . 
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III. APPLICATIONS 



A. Effect of non-additivity on the shape of depletion potentials 

1. Changes of type (a); Usa varies, but ats is fixed 

Following the procedure of we follow route (a) and demonstrate that for a given (fixed) size ratio q and packing 
fraction ■q'^'^'^ of the additive mixture the effects of non-additivity A on the shape of the depletion potential are very 
strong. 

Choosing a size ratio q = 0.1 and a fixed packing fraction in the additive mixture of rj'^'^'^ = 0.2, we can vary the 
non-additivity parameter A between A = —0.031, corresponding to the packing fraction r]l{A) of the small spheres, 
given by Eq. a, reaching the freezing packing fraction •q-^^'^'^^ — 0.494, and A — A^'-' — g, in which case ass = and 
hence 'r]l{A^'~^ = 0. Results for the depletion potentials are given in Fig. ^. For q = 0.1 and r]'^'^'^ = 0.2 the depletion 
potential obtained from DFT for additive hard spheres (A = 0) is in excellent agreement with the results of computer 
simulations |jl^ - see the comparison in Ref. [0. Moreover for A = A^^ = 0.1 we find that our calculated depletion 
potential is indistinguishable from the analytic AO result jist . 

For small degrees of positive non-additivity the main effects observed in the depletion potential are a weakening 
of the first repulsive potential barrier, due to a decreased packing fraction [given by Eq. (^)], and an increase in the 
range of attraction, i.e. the maximum of the potential shifts to larger separations h. Both effects tend to increase the 
net attraction and this is reflected in the second virial coefficient, as we shall see later. 

For negative values of A the packing fraction rj^ (A) of the small spheres increases rapidly and the contact value 
of the depletion potential increases sharply - see Fig. H. For sufficiently negative values of A the depletion potential 
can become positive at contact while the force near contact remains attractive. If the density of the small spheres 
is small enough to permit a high degree of negative non-additivity the depletion force near contact can even become 
repulsive. This is illustrated in Fig. for parameters r]^'^'^ = 0.1 and q = 0.1. Now A can take values as low as —0.060 
whilst ?7^(A) remains smaller than rf/'^^'^ . Depletion potentials for q = 0.1 and 77°'*'* — 0.3 were presented in Fig. 3 of 
Ref. PI and for q = 0.2 and rjf'^ = 0.1 were presented in Fig. 2 of Ref. ; these display similar trends with A as 
those shown here. 



2. Changes of type (b) ; aba varies, but ass is fixed 



In Fig. ^ we show the effect of changing A according to route (b). Now the packing fraction in the reservoir 77J is 
fixed at 0.2 for all values of A and ats varies. The results are very different from those in Fig. ^ which correspond to 
the same size ratio, q — 0.1, and the same 77"'''^. In the present case increasing A shifts, almost rigidly, the depletion 
potential to larger separations h leading to much deeper and longer ranged attractive wells than for A = 0. Making A 
increasingly more negative corresponds to shifting the potential to smaller separations thereby reducing the attraction 
and the height of the potential barrier. Depletion potentials for q ~ 0.2, 77J = 0.2 were presented in Fig. 3 of Ref. [ p9| ; 
these display similar trends with A as those shown here. 



B. Effect of non-additivity on the second virial coefficient 



1. Changes in B2 as a function of A 



In the subsection III A we demonstrated that introducing a rather small degree of either positive or negative non- 
additivity has a profound effect on the shape of the depletion potential. Here we investigate the effect of non-additivity 
on the second virial coefficient B2 which measures the net attraction between two big particles in the sea of small 
ones. B2 corresponds to the total effective pair potential ^effizs^r) defined in Eq. (|2l|). It follows that 



-B2 = SJ''^ + 27r / drr^{l-exp[-/3tJ2(zs;r)]} 



(36) 



with B^^ 



/3, the second virial coefficient of the pure hard-sphere system. If the depletion potential LU2{zs]r) 
generates enough attraction between the two big spheres B2 can become negative |Q . Note that B2 is a function of 
A and z,. 
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In Fig. B2 is plotted as a function of A for rj'^'^'^ = 0.2 and q = 0.1. For an intermediate value of Amin ~ 0.0279 
we find that B2 takes its minimum value of i?2(Ami„)/i3|^"^ « —0.5. The variation of B2 with A is similar to that of 
W{h = 0) = uj2{zs',crbb), the contact value of the depletion potential, although the latter has its minimum at a slightly 
lower value A « 0.016 - see inset of Fig. |7[ Provided rjl is sufhcientlty high to generate significant packing effects the 
presence of a minimum in B2 in the range > A,„i„ > A^'-' is easily understood JlSf . For a given size ratio q and 
density of small particles > the Asakura-Oosawa depletion potential given by 

W{h) = W^^'ih) = -p,d{zs)AV{h), (37) 

where Pid{zs) — p^ksT is the ideal gas pressure and AV{h) is the overlap volume excluded to the centers of the small 
spheres, h denotes the separation between the surfaces of the two big spheres so that AV{h) = for h > (Tss{0) — qo-bb- 
W^'^ is purely attractive and should always generates more net attraction than the depletion potential of an additive 
hard-sphere mixture - provided 77^ is large enough that packing effects become significant. Then the depletion 
potential consists of an attractive part close to contact and an oscillatory tail for larger separation. Packing effects 
of the small spheres reduce the range of the initial attractive part of the hard-sphere depletion potential compared to 
the Asakura-Oosawa potential and for the same value of q and we find -B^*^ < Sf'''^, at least for the parameters 
we studied. Very close to the Asakura-Oosawa limit, i.e. A < A"^*^, where the packing fraction of the small spheres, 
Eq. (j^), is small but non-zero, packing effects are minor and the depletion potential is still determined by excluded 
volume considerations. For a non-zero packing fraction, however, the pressure of the small sphere fluid is higher than 
in an ideal gas so that to first order in 77^, the virial expansion of this pressure yields 

W^^Hh) = {l+4T]^si^)W^"ih). (38) 

In Ref. |l^ we showed that this modified AO approximation is very accurate for ?7^(A) < 0.01. According to this 
approximation the depletion potential, Eq. (|3^), is more attractive than in the Asakura-Oosawa limit so that we can 
conclude that for an intermediate value of A the second virial coefficient B2 must have a minimum. It is, however, 
surprising and striking that this minimum is found to be deep, (i?|^'^ — B2{Amin)) / {B2 ^ — B2{A'^'~')) w 1.26, and 
located at a rather low degree of non-additivity. 

In the additive limit B2 is already positive for this particular mixture and a very small degree of negative non- 
additivity is sufficient to make B2 strongly positive. 



2. Changes in B2 as a function of vil and a criterion for fluid-fluid phase separation 

In an experimental situation A is not an easily controllable or tunable parameter and therefore it is most interesting 
from an experimental point of view to consider a fixed value of A and investigate the depletion potential and B2 when 
the reservoir density of the small spheres p^, a quantity that can be controlled easily in an experiment, is changed. 

If one were to take any real (asymmetric) binary mixture of hard-sphere like colloidal particles and have some means 
of determining the three interparticle pairwise potentials one could, in principle, assign three effective hard-sphere 



diameters Ubbi <^bs and Uss using standard liquid state theories |29|. In general one would not expect these diameters 
to be perfectly additive although the magnitude and sign of the non-additivity might be difficult to ascertain by any 
direct measurement. In this subsection we demonstrate that a very small degree of non-additivity reveals itself very 
clearly in the dependence of B2 on the packing fraction of the small spheres rf^. 

To this end we start with an additive mixture and plot in Fig. B2 expressed in units of the second virial coefficient 
of a pure hard-sphere system B2^ as a function of ry^ = r]'^'^'^ for various size ratios q. The qualitative behavior of 
B2 is the same for all size ratios: for small values of rf^ the reduced second virial coefficient decreases from unity in 
approximately linear fashion while for high packing fractions, rfg > 0.3, the decreasing range of attraction and the 
increasing height of the repulsive potential barrier in the depletion potential hinder a further decrease of B2 and we 
find a minimum of the second virial coefficient at roughly rf^ < 0.4 and an increase of B2 upon further increase of rf^. 

The quantitative behavior, however, depends very strongly on the value of q. For large size ratios, q > 0.2 we 
find that for all packing fractions rj^ the depletion potential cannot generate enough net attraction to make B2 
negative. This observation helps us to understand, in terms of the depletion potential, the fact that in additive binary 
hard-sphere mixtures with large size ratios q even metastable fluid-fluid phase separation does not occur . 

For smaller size ratios the minimum of B2 becomes more pronounced and B2 takes on negative values for a range 
of T]l. At a size ratio 1 : 9 (g = 0.1111) the minimum value of B2 falls below — 1.5i?|^'^. By analyzin g a large 
series of simulation results for a variety of (one-component) model fluids, Vliegenthart and Lekkerkerker [ pl| have 
shown recently that the second virial coefficient evaluated at the gas-liquid critical temperature Tc takes values which 
lie in a fairly narrow range around — 1.5-B|^'^. For the model fiuids considered in Ref. pl| gas-liquid coexistence 
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can only occur if 62/62^ < — 1.5. If we assume that this empirical criterion is applicable to the effective one- 
component system described by our pairwise potentials ^eff it follows that only systems with 32/82^ lying below 
the horizontal line in Fig. |^ could exhibit (metastable) fluid-fluid coexistence. It is important to emphasize that the 
criterion is empirical and that it was developed for model fluids in which the attractive part of the pairwise potential is 
monotonically increasing with interparticle separation r, unlike our present effective potentials. Moreover the criterion 
does not predict whether the gas-liquid coexistence is stable or metastable w.r.t. fluid-solid coexistence. Recall that 
the shorter range of the attractive potential the more likely is the gas- liquid transition to become metastable 
In a simulation study of the effective one-component Hamiltonian for an additive hard-sphere mixture, metastable 
fluid- fluid phase separation was found for g = 0.1 and q = 0.05 Q. For q = 0.1 fluid- fluid coexistence occurred for 
rjl > 0.29, whereas for q = 0.05 this occurred for ry^ > 0.165. These results are in keeping with the predictions of 
the empirical criterion. Note that the second intersection of 62/62^ with the horizontal line for g = 0.1 in Fig. ^ 
suggests a possible upper critical point near -ql = 0.43 |29|. 

If we introduce a very small degree of positive non-additivity according to route (a), i.e. we keep atb and abs 
constant so that crss(A) becomes smaller than for the additive case (see Eq. we find a dramatically different 
situation. In order to be of relevance to an experimental situation, where it might be impossible to rule out a small 
degree of non-additivity, we set A = g/20 for each choice of q. 

As expected, positive non-additivity leads to a slight increase in the width of the depletion layer and therefore more 
net attraction than in the additive case. The behavior of B2 as a function of 77^ (A), however, changes qualitatively as 
can be seen in Fig. ||. For all size ratios considered here the second virial coefficient becomes negative and 62/62^ 
falls below the —1.5 line prior to freezing of the small particles. The smaller q, the smaller the value of r]l when this 
line is crossed. There is no minimum in B2. Thus according to the empirical criterion all the mixtures considered 
here should exhibit (metastable) fluid-fluid coexistence. Once again we cannot say whether this transition is stable 
w.r.t. the fluid-solid transition. 

Our results provide some understanding, in terms of the depletion potential description, of why only small degrees 
of positive non-additivity might lead to fluid-fluid phase separation in asymmetric binary hard-sphere mixtures. We 
recall that Biben and Hansen |^ found for g = 0.1, on the basis of the Barboy-Gelbart |^ equation of state, that a 
value A = 0.01 was sufficient to produce a fluid-fluid transition at a total packing fraction < 0.5. Later Dijkstra |Tc| ] 
carried out a series of Gibbs ensemble Monte-Carlo simulations of the binary mixture for q = 0.1 and varying degrees 
of positive non-additivity. She found that it was possible to have fluid-fluid demixing for a total packing fraction 
< 0.5, provided A was sufficiently large. 

On the other hand, introducing a small degree of negative non-additivity [again via route (a)] decreases the width of 
the depletion layer compared to the additive case so that the net attraction should also decrease. We set A = —q/20 
for each q and find that these small negative values of A are sufficient to change significantly the shape of B2 versus 
?7^(A) from that of the additive case. 

This is illustrated in Fig. |l^ where we show that for size ratios q > 0.2 the second virial coefficient changes little 
as function of 77^ (A). For large values of rjl we even find that B2 > 82^. For a size ratio of q ~ 0.1, which in the 
additive case was sufficiently small for B2 to become strongly negative, B2 remains positive for all packing fractions. 
Only for very small q can the depletion potential in non-additive mixtures with negative A generate sufficient net 
attraction to drive B2 negative. 

We conclude this discussion by emphasizing how sensitive B2 is to changes in the depletion potential. In Fig. |ll| 
we show the depletion potentials calculated for q — 0.1 and fixed ?7s(A) = 0.3 for three values of A. Although the 
three potentials appear rather similar they yield very different values of B2. From Figs. p|-p^ we see that B2I B2^ 
is about —1.21 for A = 0, is strongly negative {B2/ B2^ < —9) for A = 0.005 (the empirical criterion would predict 
fiuid-fiuid phase separation) and is positive {B2/B2^ — 0.36) for A = —0.005. 



IV. DISCUSSION 



The main results of our study of the equilibrium statistical mechanics of non-additive binary hard-sphere mixtures 
may be summarized as follows: 

1. The formal technique of integrating out the degrees of freedom of the small spheres in order to obtain an effective 
Hamiltonian H'^^-^ for the big spheres is equally valid for non- additive as for additive mixtures, where it has proved 
particularly useful |^ for determining the phase behavior of asymmetric systems. We have provided expressions for 
zero, ilo, one, fli, and two body contributions to H'^^f which can be evaluated in simulations of the small sphere fluid 
(Sec.^LAl). 

2. We showed that the same contributions are readily calcu lated using the fundamental measure DFT of Rosenfeld, 
a theory developed originally for additive hard-spheres (Sec. [IB). By calculating Oq and fii we were able to derive 
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an explicit approximation ( p4[ ) for the packing fraction rjs of small spheres in the mixture at given reservoir fraction 
r]^ and big-sphere fraction rjb- Comparison with a previous simulation study for additive mixtures with small size 
ratios q showed that the approximation is very accurate (Fig. ||) . 

3. The two body contribution to H'^^^ is a sum of effective pairwisc potentials ^eff{zs; r) — <i>bb(r) +lu2{zs', f)- Our 
DFT approach provides a powerful means of determining the depletion potential lo2{zs] r) for non-additive mixtures. 
Provided three and higher body terms are small, as is expected for small q, it is solely uj2 that determines big-big 
correlation functions and the phase equilibria of the mixture. 

4. We described two different routes to introducing non-additivity A and gave examples of the depletion potentials 
and the second virial coefficient B2 associated with the corresponding effective potential ^e// obtained from both 
routes (Sec. HI). What emerges is that the depletion potential depends strongly on whether route (a)) a^s fixed, ct^s 
varies or route (b) cjgg fixed, Cfts varies, is employed. B2 is a sensitive indicator of the shape and range of the effective 
potential and exhibits considerable variation with q, A and t]^. (Figs |7|-[lC|) . 

5. On the basis of the empirical criterion B2/B2^ < —1.5 [pl|, we showed that fluid-fluid phase separation is much 
more likely to occur for a small degree of positive non-additivity, A > 0, than in additive mixtures, A = 0, with 
the same size ratio. Our results provide a guide to which binary hard-sphere mixtures might exhibit a stable, w.r.t. 
fluid-solid, fluid-fluid separation and we hope that these might stimulate further computer simulation studies. 

We conclude by turning to the question of whe ther the s trong effect of non-additivity on the depletion potentials and 
on the virial coefhcients found in sections [II A and III B has any implications for experiments on colloidal systems. 

One-component colloidal suspensions which mimic very closely a hard-sphere system can be created because any 
small residual short-range interactions remaining after refractive index matching are very well approximated by hard 
spheres with an effective hard-sphere diameter p3] . But creating a truly additive binary colloidal suspension is much 
more difficult, since this implies an additional constraint on the value of the effective hard-sphere diameters, namely 
that 2ahs — o'bt + ctss- The small residual interactions in an experimental system designed to mimic binary hard- 
sphere mixtures can easily introduce non-additivity |29f|; non-additivity is probably the rule and perfect additivity 
the exception. For example, in an earlier publication |[L1[, one of us has shown by a simple argument that for both 
sterically and electrostatically stabilized asymmetric binary colloids, 2abs is likely to be smaller than (crht + o'ss)-, which 
implies a small negative non-additivity. This in turn suggests that the well-depth at contact W{h = 0) is smaller 
than what would be expected for an additive system. For short-ranged depletion systems the location of a fluid-solid 
liquidus line can be roughly correlated to W{h = 0) ]lT| , p2[ ; one would therefore expect the experimental liquidus line 
to occur at larger values of 77^ than what is predicted for a purely additive binary hard-sphere system. Experimental 
results do seem to follow this trend Q. However, since the experimental phase boundaries are typically plotted with 
the small-particle packing fraction 77^ on the y-axis, it is not always clear whether discrepancies with the additive 
theory arise from non-additivity, or from small errors in the measurement of ass, which enters rjs as a^^. 

Instead of focusing on phase boundaries, we propose that direct measurements of the osmotic second virial coefficient 
B2 as a function of rjs should be a much more sensitive measure of the existence of non-additivity, and may even 
provide an independent way to ascertain the value of abs, which is otherwise very hard to determine. As illustrated 
in Figs H, ^ and |l^, different degrees of non-additivity induce clear qualitative differences in the dependence of B2 
on r/s, implying that one does not require a high level of quantitative accuracy in measurements of B2 in order to 
distinguish clearly between negative and positive non-additivity. 

For colloidal suspensions B2 is typically extracted from the low density limit of the osmotic equation of state, 
measured by static light scattering. This is non-trivial since it requires extracting the contribution from the big 
particles to the total scattering intensity. Such measurements were first carried out by de Hek and Vrij in 1982 |3|] 
for a colloid-polymer mixture with a size ratio 2Rg/acc ~ 1 (here Rg is the radius of gyration of the polymers); they 
found a clear trend towards negative values of B2 upon increasing the polymer concentration, which is con sisten t with 
the expected positive non-additivity in such colloid-polymer systems. In contrast, the results in section [II B imply 
that for an additive binary colloid mixture, size-ratios of ass/o^bb ^0.2 are required to drive B2 negative. For negative 
non-additivity even smaller size-ratios are required. Such qualitative effects should be visible in experiments. 

We note in passing that the effect of increasing the polymer concentration on the second osmotic virial coefficient of 
a globular protein-polymer solution has been measured recently | |36[ |. However, these experiments are in the protein- 
limit 2Rg/acc S> 1, where the concepts of negative and positive non-additivity are less useful. 
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APPENDIX: LOW DENSITY LIMIT OF i^i 



Here we demonstrate that the one-body term Puji{zs), given in Eq. reduces in the Umit ^ to the exact 
low density Umit, i.e. 

lim I3ui{z,) =^ p^Vb+s, (Al) 

with Vb+s = 47rtT^^/3 the volume of a spherical cavity of radius abs which is excluded to the centers of the small 
spheres. Recall that 0Ji{zs) is the excess chemical potential of a big hard sphere in a sea of small ones. 

In order to take the low density limit we note that to leading order in powers of pi the equation of state, Eq. ([2^), 
reduces to 



the surface tension, Eq. (p2|), to 



]imf3p = p:, (A2) 



lim/37=^p:, (A3) 



the coefhcient of the term in Plui linear in Rb reduces to 

hm , 



and, finally, 



^lim^-ln(l-r,:) = -a,^,p:. (A5) 



It follows that 

47r „o ^ . ass ^ . „ / cr.,c \ 2 ^ 47r / fjc, \ 3 



\hn^Pu;, = ^RIp: + ^.Rf-fpl + 4nRb (^) p^ + ^ (^) p^. (A6) 
Recalling that Rb = abs — c^ss/^ we find limp^^o P^^i ~ ^'^tsP^s which is Eq. ( jAll ). 



[1] T. Biben and J.-P. Hansen, Phys. Rev. Lett. 66, 2215 (1991). 

[2] M. Dijkstra, R. Van Roij, and R. Evans, Phys. Rev. Lett. 81, 2268 (1998); ihid 82, 117 (1999); Phys. Rev. E 59, 5744 
(1999). 

[3] S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 (1954). 

[4] S. Asakura and F. Oosawa, J. Pol. Sci. 33, 183 (1958); A. Vrij, Pure and Appl. Chem. 48, 471 (1976). 
[5] A. P. Cast, C. K. Hall, and W. B. Russel, J. Coll. Int. Sci. 96, 251 (1983). 

[6] H. N. W. Lekkerkerker, W. C. K. Poon, P. N. Pusey A. Stroobants, and P. B. Warren, Europhys. Lett. 20, 559 (1992). 
[7] M. Dijkstra, J.M. Brader, and R. Evans, J. Phys.: Cond. Matter 11, 10079 (1999). 
[8] E.J. Meijer and D. Frenkel, J. Chem. Phys. 100, 6873 (1994). 
[9] T. Biben and J.-P. Hansen, Physica A 235, 142 (1997). 
[10] M. Dijkstra Phys. Rev. E 58, 7523 (1998). 

[11] A. A. Louis, R. Finken, and J.-P. Hansen, Phys. Rev. E, 61, R1028 (2000). 

[12] In Ref. [0 a truncated virial expansion of the Derjaguin approximation to the depletion potential for additive hard-sphere 
mixtures (see B. Gotzelmann, R. Evans, and S. Dietrich, Phys. Rev. E 57, 6785 (1998)) is employed with an effective size 
ratio q^ff = (?(1 -|- A) -|- A to take account of non-additivity. 

[13] T. Biben, P. Bladon, and D. Frenkel, J. Phys.: Condens. Matter 8, 10799 (1996). 

[14] R. Dickman, P. Attard, and V. Simonian, J. Chem. Phys. 107, 205 (1997). 

[15] See, e.g., J.C. Crocker, J.A. Matteo, A.D. Dinsmore, and A.G. Yodh, Phys. Rev. Lett. 82, 4352 (1999); D. Rudhardt, C. 
Bechinger, and P. Leiderer, Phys. Rev. Lett. 81, 1330 (1998); and references therein. 



13 



R. Roth, R. Evans, and S. Dietrich, Phys. Rev. E 62, 5360 (2000). 

B. Gotzelmann, R. Roth, S. Dietrich, M. Dijkstra, and R. Evans, Europhys. Lett. 47, 398 (1999); R. Roth, doctoral thesis, 

Bergische Universitat Wuppertal (1999). 

R. Roth and R. Evans, Europhys. Lett. 53, 271 (2001). 

J. R. Henderson, Mol. Phys. 50, 741 (1983). 

M. Dijkstra, R. van Roij, and R. Evans, J. Chem. Phys. 113, 4799 (2000). 

T. L. Hill, in An Introduction to Statistical Thermodynamics, Dover publications. Inc., New York (1986). 
A. A. Louis, Phil. Trans. Roy. Soc. A 359, 939 (2001). 

D. J. Goulding, Ph.D. Thesis University Cambridge (2000); D. Goulding and S. Melchionna, to appear in Mol. Phys.. 

S. Melchionna and J.-P. Hansen, Phys. Chem. Chem. Phys. 2, 3465 (2000). 

Y. Rosenfeld, Phys. Rev. Lett. 63, 980 (1989); J. Chem. Phys. 98, 8126 (1993). 

R. Evans, Adv. Phys. 28, 143 (1979). 

Y. Rosenfeld, M. Schmidt, H. Lowen, and P. Tarazona, J. Phys.: Condens. Matter 8, L577 (1996); Phys. Rev. E 55, 4245 
(1997). 

H. N. W. Lekkerkerker and A. Str oobants, Physica A 195, 387 (1993). 



A. A. Louis and R. Roth, preprint |cond-mat/0102049| (2001). 

As first pointed out by H. De Hek and A. Vrij, J. Coll. Int. Sci. 84, 409 (1981). 

G. A. Vliegenthart and H. N. W. Lekkerkerker, J. Chem. Phys. 112, 5364 (2000); see also M. G. Noro and D. Frenkel J. 
Chem. Phys. 113, 2941 (2000). 

B. Barboy and W. M. Gelbart, J. Chem. Phys. 71, 3053 (1979). 

P. N. Pusey and W. van Megen, Nature (London) 320, 340 (1986). 

J.S. van Duijneveldt, A.W. Heinen, and H.N.W. Lekkerkerker, Europhys. Lett. 21, 369 (1993); A. D. Dinsmore, A. G. 
Yodh, and D. J. Pine, Phys. Rev. E 52, 4045 (1995); A. Imhof and J. K. G. Dhont, Phys. Rev. Lett. 75, 1662 (1995). 

H. De Hek and A. Vrij, J. Coll. Int. Sci. 88, 258 (1982). 

A. M. Kulkarni, A. P. Chatterjee, K. S. Schweizer, and C. F. Zukoski, Phys. Rev. Lett. 83, 4554 (1999). 



14 




FIG. 1. The centers of the small spheres are excluded from the shaded depletion layer around each big sphere of diameter 
abb- When changing the non-additivity A by route (a), the big-small diameter abs is kept constant. For a fixed abb this implies 
that the depletion layer thickness I = abs — crbb/2 is also held constant. The small sphere diameter crss(A) is decreased for 
A > and increased for A < 0, and since the number density pi is fixed, the small sphere packing fraction r]s{A) also changes. 
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A>0 A<0 




FIG. 2. When changing A by route (b), the big-small diameter diameter ats changes. For a fixed att this implies that the 
depletion layer thickness / = abs — crbb/2 also changes. The small-particle diameter ass is kept constant, and since is fixed, 
the small sphere packing fraction rjs is also constant. 
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FIG. 3. The small sphere packing fraction rjs of an additive binary hard-sphere mixture, with size ratio (a) 0.10 and (b) 
0.05, versus that of the reservoir rjl for several big sphere packing fractions rjt. The squares denote the direct simulation data of 
while the lines denote the results of Eq. (|^. Note the significant deviations from linearity in both theory and simulations. 
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h / 033(0) 

FIG. 4. The depletion potential W{h) = uJ2[zs;(Jbb + h) between two big hard spheres in a sea of small hard-spheres 
calculated for a size ratio q = Q.l and a range of non-additivities A treated according to route (a). The number density of 
the small spheres, pj, is fixed with packing fraction 77"'^'' = pl'^{(Tsf'^)^ /& = 0.2. A = corresponds to an additive hard-sphere 
mixture, h is the separation between the surfaces of the big spheres and ass{Q) = qcrhb- 
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h / 033(0) 

FIG. 5. As in Fig. ^ but now the packing fraction 77°'*'^ is fixed at 0.1. For these negative values of the non-additivity A the 
depletion potential is repulsive near contact and for A = —0.05714 the depletion force is repulsive near contact. 
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FIG. 6. The depletion potential W{}i) = ui2{zs;o-bb + h) between two big hard spheres in a sea of small hard-spheres 
calculated for a size ratio q = 0.1 and a range of non-additivities A treated according to route (b). The packing fraction in 
the reservoir rjl = 0.2 remains constant for all values of A. A = corresponds to an additive hard-sphere mixture, h is the 
separation between the surfaces of the big spheres and (Tss(O) = qatb- These results should be contrasted with those in Fig. ^ 
- note the difference between the vertical scales. 
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FIG. 7. The second virial coefficient of the big spheres B2 in units of = 2Tvalf,/3 as a function of the non-additivity A. 
These results correspond to the depletion potentials of Fig. ^ i.e. r;"'*'^ = 0.2 and q = 0.1, for which the maximum non-additivity 
is A'*'' = 0.1. The inset shows the contact value l3W{h = 0) = (3u)2{zs\ Ubb) of the depletion potential as a function of A. 
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FIG. 8. The reduced second virial coefficient of the big spheres in additive hard-sphere mixtures (A = 0) for various size 
ratios q versus •q^'^'^ , the packing fraction of smaU spheres in the reservoir. According to the criterion of Ref. |3l|, (metastable) 
fluid-fluid phase separation can only occur if B2/B§^ < —1.5 (horizontal line). 
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FIG. 9. The reduced second virial coefficient of tlie big spheres in hard-sphere mixtures with a small positive non-additivity 
A — +q/20 for various size ratios q versus ril{A), the packing fraction of small spheres in the reservoir. In contrast to the 
additive case, Fig. ^, 82/82^ always falls below —1.5 indicating that (metastable) fluid-fluid phase separation could occur for 
all size ratios shown here. 
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FIG. 10. The reduced second virial coefficient of the big spheres in hard-sphere mixtures with a small negative non-additivity 
A = — g/20 for various size ratios q versus ??J(A), the packing fraction of small spheres in the reservoir. Note that B^jB^^ 
falls below —1.5 only for the smallest ratio, g = 0.05, considered here. 
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FIG. 11. Depletion potentials for q = 0.1 and fixed t^KA) = 0.3 for A = -0.005, and +0.005. Althougli these potentials 
do not differ drastically from each other the corresponding B2, see Figs. p[-px|, take very different values: = 0.36, 

-1.21 and < -9 for A = -0.005, 0, and +0.005, respectively. 
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